Correlations in avalanche critical points 
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Avalanche dynamics and related power law statistics are ubiquitous in nature, arising in phe- 
nomena like earthquakes, forest fires and solar flares. Very interestingly, an analogous behavior is 
associated with many condensed matter systems, like ferromagnets and martensites. Bearing it in 
mind, we study the prototypical 3D RFIM at T = 0. We hnd a Hnite correlation between waiting 
intervals between avalanches and the previous avalanche size. This correlation is not found in other 
models for avalanches, such as the standard BTW model, but it is experimentally found in earth- 
quakes and in forest fires. Our study suggests that this effect occurs in critical points which are at 
the end of an athermal first-order transition line separating two behaviors: one with high activity 
from another with low activity. 

PACS numbers: 05.70.Jk, 05.40.-a, 75.60.Ej,, 75.40.Mg, 75.50.Lk 



In the last few years much experimental and theoreti- 
cal effort has been devoted to the study of avalanche pro- 
cesses. Deep understanding of the statistical correlations 
in such stochastic processes is needed in order to make 
advances towards predictability. The importance of the 
subject is beyond discussion due to the many implications 
in natural disasters and social crises. Avalanche pro- 
cesses are characterized by extremely fast events whose 
occurrences are separated by waiting intervals without 
activity. The magnitudes characterizing avalanches (en- 
ergy, size, duration) are, in most cases, statistically dis- 
tributed according to a power law p{s)ds s~'^ds char- 
acterized by a critical exponent t: extremely large events 
hardly occur, whereas small events are very common. 
This is the famous Gutenberg-Richter law for earth- 
quakes. Power-law distributions have been found not 
only for other large-scale natural phenomena ranging 
from solar flares [l[ to forest-fires 0], but also in labora- 
tories associated with many condensed matter systems: 
condensation Q, ferromagnets martensitic transi- 
tions [1], superconductivity etc. However, statis- 
tics of waiting intervals 6 are not often studied, espe- 
cially in condensed matter. The distribution p{S)dd has 
been described by different laws including exponentials 
and also power-laws. For solar flare statistics and earth- 
quake models, it has been found 0, Hi that the unavoid- 
able threshold definition (separating activity from inac- 
tivity) alters the distribution p{5) and that this thresh- 
olding effect is a signature of the existence of correla- 
tions 0, [13] ■ Direct measurement of the correlations 
between waiting intervals and avalanche sizes has been 
obtained by measuring the constrained interval distribu- 
tions Pprev{S\s > So) and Pnext{S\s' > So)- Thesc are the 
probabilities of having a waiting interval 6, given that the 
previous (s) or the next (s') avalanche is larger than sq. 
For earthquake and forest-fire statistics, while Pnext has 
been found to be independent of Sn, Ppr-e.v does exhibit 
significant changes when varying sq 



In this letter we study some statistical correlations 
for the 3D-RFIM at T = with metastable dynamics 
based on the local relaxation of single spins. This model 
was introduced Il3| for the study of Barkhausen noise 
in ferromagnets [J| and acoustic emission in martensitic 
transitions [1], and has been used as a prototype for the 
study of crackling noise and other avalanche phenomena 
l4 , ll5|. The model is defined on a cubic lattice with 
size N = L^. At each lattice site there is a spin variable 
Si = ±1 {i — 1, . . . ,N) that interacts with its nearest 
neighbors (n.n.) according to the Hamiltonian (magnetic 
enthalpy): 
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The local random fields hi are Gaussian distributed with 
zero mean and standard deviation a. This parameter not 
only allows the critical behavior [a — ac — 2.21 ±0.01) to 
be studied but also subcritical (cr > Uc), and supercritical 



(cr < (Jc) regimes [16|, |17[. This tuning is absent in so- 



called Self-Organized Criticality models like the original 
BTW model [18|, in which the critical state is reached 
after waiting for a certain time without any parameter 
adjustment. In the RFIM, the time variable is replaced 
by the external field H , which is adiabatically increased 
from — oo to oo. The system responds by an increase of 
the order parameter m = X]i=i Si/N (magnetization per 
spin) from —1 to 1. The spins flip according to the dy- 
namical rule Si — sign(^,- Sj + hi + H) (the first sum 
runs over the n.n. of spin Si), which corresponds to a 
minimization of the local energy. This non-equilibrium 
dynamics leads to hysteresis and avalanches when many 
spins fiip at constant field. The avalanche size s is de- 
fined as the number of spins flipped until a new stable 
state is reached. Avalanches are separated by field wait- 
ing intervals S without activity. When a — ac and H is 
close to He — 1.43 ± 0.05 [H, [3], the avalanche size dis- 
tribution becomes approximately a power lawp(s) ^ s""^ 



2 



(with r ~ 1.6 [l^). The properties of large (percolating) 
avalanches at the critical point have been previously dis- 
cussed [l^. Critical avalanches have been found to be 
fractal so that (s)c ~ L'^-f with df ~ 2.88. When a > ac, 
the avalanche distribution is exponentially damped, i.e. 
all avalanches are negligible compared to system size, and 
the magnetization m evolves continuously when the sys- 
tem is infinite. For a < avalanches are also infinitesi- 
mally small compared to L^, except for a unique infinitely 
large and compact (s ^ L^) avalanche corresponding 
to a first-order phase transition between a phase with 
low m and a phase with high m. The first-order tran- 
sition line can be linearly approximated by the equation 
Htia) =Hc[l- B'{a - ac)/ac] with B' = 0.25. [la] 

In this letter we will concentrate on the analysis of the 
field waiting intervals S, and their correlation with the 
previous and next avalanche sizes. We have first checked 
that intervals are exponentially distributed. Fig. [T] shows 
the distribution of waiting intervals for a = 2.21, L = 30 
and three different if-ranges. One can observe that be- 
fore and after the transition region the distributions are 
very well described by exponentials (continuous lines), 
whereas in the transition region the distribution becomes 
a linear mixture of exponentials. The mean value (6), 
which is the only parameter characterizing such distribu- 
tions, depends on a and H, and below ac it exhibits 
a discontinuity AS when H increases and crosses the 
first-order transition line from the region of high activity 
(small (S)) to the region of low activity (large {6) ), as 
shown in the inset of Fig. [TJ We can, therefore, compare 
the behavior of the discontinuity AS to that of an order 
parameter. Note that, given the finite size of the system, 
the pseudo-critical point where the mixture of exponen- 
tials becomes a single exponential distribution will be 
located at ct > 2.21. For this reason, although data in 
Fig. [T] correspond to dc, one can still observe a range of 
fields with the two-peak distribution. 

Besides the histogram analysis, we have numerically 
checked that for all a and H, (S) - VWT^^W.H 
which is more evidence of the exponential character of 
p{5). The continuous line in the inset of Fig. [1] shows 
this agreement, except for some deviations close to the 
transition line due to finite-size effects. 

For the following discussions it is interesting to ana- 
lyze the finite-size dependence of (S). A simple argu- 
ment can be used to state that far from the critical point, 
since the correlation length is finite, the probability for 
an avalanche to start when the field is increased by dH is 
proportional to the number of triggering sites and thus to 
L^. This implies that (S) ~ L~^. Fig. [2] shows examples 
of this behavior for different values of a and H. At crit- 
icality, the correlation length diverges and the triggering 
argument may be too naive. Nevertheless, as shown in 
Fig. [21 numerical simulations indicate that the exponent 
is always very close to 3. Uncertainties in He and <Tc do 
not allow for an accurate enough finite-size scaling anal- 
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FIG. 1: Log-log plot of the histograms corresponding to the 
distribution of waiting intervals S for L = 30, a — 2.21 and 
three different field ranges, indicated by the legend. Note that 
the horizontal scale is in dB and bins are logarithmic. The 
dashed lines show fits corresponding to the exact exponential 
behavior. Data have been obtained by averaging over 50000 
realizations of disorder. The inset shows the behavior of (S) 
(symbols) and ^/{S^)~-~{S)^ (continuous line) as a function of 
H foT a = 1.95 and L = 30. 



ysis to determine small variations. For the discussions 
here the exact value of this exponent is not needed and 
we will assume that (5) with z ^ 3. 
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FIG. 2: Log-log plot of the average waiting interval (S) as a 
function of the system size L for different values of a and H 
as indicated by the legend. The continuous line shows the 
behavior L~'^ and the discontinuous lines are fits to the last 
2-3 points for each series of data. Data are averaged over 
50000 disorder realizations. The inset shows the behavior of 
(s^) — for the values of a and H indicated. The contin- 
uous lines show the behavior L'^f (above) and L^^^ (below). 



Let us now focus on the study of correlations and con- 
sider a sequence of two avalanches, the first with size 
s, then a waiting interval S, and the next avalanche with 
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size s' . We define the following two correlation functions: 

{sS)-{s){S) 



Ps,S 
PS,s' 



(2) 
(3) 



Fig. [3] shows examples of the behavior of the correlation 
functions for different system sizes and a — ac — 2.21 as a 
function of the scaling variable v = [{H — He) /He] L^^^ 
that measures the distance to the critical point. The 
exponent /x = 1.5 has been taken from the literature [l^. 
(Note that this scaling variable only takes three values in 
the thermodynamic limit v = ±oo, 0). 

For the understanding of the behavior of the correla- 
tion functions for increasing system size, one must first 
discuss what the expected behavior of a generic correla- 
tion function p{H, a = CTc, L) close to H,, is likely to be. 
It should be borne in mind that correlations are, by defi- 
nition, bounded between —1 and 1. Therefore no critical 
divergences can occur with increasing L. Consequently, 
any critical correlation either goes to zero or tends to 
a constant value. In this second case, it should exhibit 
scaling behavior p{H, CTc, L) ~ p{(Jc, v). 

The first observation from Fig. [3] is that ps^s' is much 
smaller than p^^s not only at the critical point [v = 0), 
but also for other values of the field. In addition, the 
data shows that ps^s' systematically decreases in abso- 
lute value with increasing system size. Therefore the nu- 
merical data is consistent with a vanishing ps,s' in the 
thermodynamic limit. 

The second important observation is that correlation 
between an avalanche size and the next waiting time, ps^s 
exhibits a constant value ~ 0.4 at the critical point. The 
overlap of the curves is rather good, specially if one takes 
into account the fact that even at the critical point there 
are non-critical avalanches which may slightly perturb 
the scaling function behavior. The fact that the scaling 
function in Fig. [3{a) goes to zero for v ±oo indicates 
that the finite correlation only survives exactly at the 
critical point, but vanishes for fields both above and be- 
low. We would like to note that the result of a finite cor- 
relation ps^s is not in contradiction with the Poissonian 
character for the triggering instants of the avalanches. 

A second way to go deeper into the understanding of 
the nature of the correlations is the direct measurement 
of the restricted distribution of intervals Pprev{S\s > So) 
and Pnext{5\s' > So). Fig.|3]shows, as an example, the de- 
pendence of these two distributions for L = 30, cr = 2.21 
and 1.40 < H < 1.45. The distribution of intervals 
Pprev{S\s > So) clcarly exhibits a dependence on the pre- 
vious avalanche size, whereas PnextiS\s' > So) is indepen- 
dent of So- Thus, the larger the size of an avalanche, the 
larger the probability that the following waiting inter- 
val is large. We should note at this point that a similar 
same causal dependence has been found for the statis- 
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FIG. 3: Correlation functions defined by Eq. ((2)1 and (|3} as 
functions of the scaling variable for g = 2.21 and different 
system sizes as indicated by the legend. Data correspond to 
averages over 50000 disorder realizations and field intervals 
with size Aiif = 0.005. Note the different vertical scale in the 
two figures. 
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FIG. 4: Log-log plot of the restricted distribution of intervals 
5 given that the previous (next) avalanche is larger than so. 
Ppreu has been displaced two decades up in order to clarify 
the picture. 



11, nj, but with a 



tics of earthquakes and forest fires 
different sign. The larger the size of an earthquake, the 
smaller the waiting time to the next event. 

The origin of such correlations in our RFIM can be 
understood by noting that, for a finite system below ctc, 
after an avalanche with a large size (of the order L'^), 
one can "guess" that the system has jumped to the low- 
activity region and, therefore, one can "predict" that the 
next interval 5 will be large. We can provide an heuris- 
tic argument of why such finite-size correlations vanish 
in the thermodynamic limit, everywhere except at the 
critical point. Let us analyze the behavior of the nu- 
merator and the two square roots in the denominator 
in Eq. [21 First note that, given the fact that p{5) is 
exponential, the fluctuations of (5^) — (5)^ behave as 
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6 ^ L~^. The fluctuations of s behave differently below 
CTc and at CTc. If we consider an interval AH that crosses 
the first-order transition region, below ac we will have 
a number of avalanches proportional to contributing 
to the averages. Among these, most will display a small 
size ~ but one will have a size ~ L^. When comput- 
ing (s) we get a (s) ^ behavior, but when computing 
\/ (s^) — (s}2 we will get L'^/^, as shown in the inset of 
Fig. [2 for a = 1.95 and the interval 1.465 < H < 1.470 
crossing the transition line. Exactly at criticality, the dis- 
tribution of avalanches becomes a power-law. This means 
that the averages (s) and (s^) should be computed by in- 
tegrating the distribution from 1 to the largest avalanche 



max 
f 



and (s^) Li 



By 



2dt 



size which has a fractal dimension s 
integration, one trivially gets (s) ^ L 
The fluctuations therefore will go as ^ (s^) — (s)^ ~ L"'! , 
as shown in the inset of Fig. [2l 

To study the numerator in Eq. ([2]) one should note that 
{s5)-{s){8) = {s'5')-{s'){5) = {s'{5'^5)). This average 
measures the avalanches that carry an associated change 
in 5. For ct < ctc, among the set of avalanches in a 
interval Ai7, there is one avalanche with size Li^ associ- 
ated with a change A(5 ^ . The rest of the avalanches 
have a size L*^ and carry no change in 5. Therefore the 
numerator in Eq. [2] goes as L^^ and consequently, below 
CT„ p{s,5) - L-y{L^/^L-'-) - L-3/2 ^ 0. 

For a — ac the behavior of the numerator in ^ is much 
more intricate since different kinds of critical avalanches 
exist close to the critical point. Apart from a number 
of non-critical avalanches, there is an infinite number 
L^) of spanning avalanches and another infinite num- 
ber ~ ij^rzBa Qf critical non-spanning avalanches. It is 
difficult to argue which of these avalanches have an asso- 
ciated AS. We cannot provide a definite argument, but, 
at least, the product {s){6) (appearing in the numerator 
of Eql2]) behaves as L'^^ L~^. Therefore it is plausible that 
at (Tc, p{s,S) ~ L''f-''/{L'^fL-^) - 1, which justifies the 
finite value of the correlation p{s, S) found numerically, 
as shown in Fig. [3l 

In order to compare with the avalanche models based 
on SOC, we have performed an analysis of the same prob- 
ability densities and correlations for the 2D-BTW model. 
In this case, the waiting times are discrete since they are 
identified with the number of grains added before a new 
avalanche starts. For large enough systems, these wait- 
ing intervals are distributed according to the geometric 
distribution which is the discrete version of the exponen- 
tial distribution and all the two correlation functions ^ 
and ^ clearly vanish. 

In summary, we have numerically studied the T ~ 
RFIM with metastable dynamics as a prototype for 
avalanche processes in condensed matter systems which 
display an underlying first-order phase transition. The 
sequence of avalanches and waiting times can be con- 
sidered as a compound Poisson process since waiting in- 



tervals tend to be exponentially distributed. Neverthe- 
less, correlations between the avalanche size and the next 
waiting time exist at the critical point in the thermo- 
dynamic limit. Such a causal correlation p{s,S) ^ 
has been found experimentally in earthquakes and forest 
fires, although with a different sign. This difference could 
easilly be explained by considering a distonintuity AS < 
0. An experimental challenge for the future is to look 
for these effects in laboratory experiments on condensed 
matter systems exhibiting avalanches (Barkhausen noise, 
acoustic emission in martensites, etc). 
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